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(-H The difficulty of simulating quantum systems, well-known to quantum chemists, prompted the 

O , idea of quantum computation. One can avoid the steep scaling associated with the exact simula- 

1^ tion of increasingly large quantum systems on conventional computers, by mapping the quantum 
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I on small quantum information processors. 
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1 INTRODUCTION 

One of the greatest challenges in quantum chemistry is to fully understand the 
complicated electronic structure of atoms and molecules. Over the last century, 
enormous progress has been made in describing the general behavior of relatively 
simple systems. In particular, combined with physical insights, elegant computa- 
tional approaches, ranging from wavefunction methods to quantum Monte Carlo 
and density functional theory, have been developed. The challenge is that the 
Hilbert spaces of quantum systems grow exponentially with system size. There- 
fore, as these methods are extended to higher accuracy or to larger systems, the 
computational requirements become unreachable with current computers. This 
problem is not merely a consequence of technological limitations, but stems from 
the inherent difficulty of simulating quantum systems with computers based on 
classical mechanics. It is therefore important to know if the computational bottle- 
necks of classical computers can be solved by a computing model based on quan- 
tum mechanics — quantum computation — whose development has revolutionized 
our understanding of the connections between computer science and physics. 
The idea of mapping the dynamics of a quantum system of interest onto the dy- 



namics of a controllable quantum system was proposed in 1982 by Feynman (43) 



and developed in 1996 by Lloyd (84). Such a quantum computer would be able 
to obtain information inaccessible with classical computers. Consequently, quan- 
tum simulation promises to be a powerful new tool for quantum chemistry. In 
this article, we review the recent applications of quantum simulation to chemical 
problems that have proven difficult on conventional computers. After introduc- 
ing basic concepts in quantum computation, we describe quantum algorithms 
for the exact, non-adiabatic simulation of chemical dynamics as well as for the 
full-configuration-interaction treatment of electronic structure. We also discuss 
solving chemical optimization problems, such as lattice folding, using adiabatic 
quantum computation. Finally, we describe recent experimental implementations 
of these algorithms, including the first quantum simulations of chemical systems. 

2 QUANTUM COMPUTATION 

2.1 Differences between quantum and classical computation 

There are fundamental differences between quantum and classical computers. 
Unlike the classical bit, which is always either a '0' or a '1', the basic unit 
of quantum information is the qubit (Fig. [I]), which can be in a superposi- 
tion of |0) and |1): a|0) -|- f3\l). States of n qubits are elements of an ex- 
ponentially large, 2^-dimensional, Hilbert space, spanned by a basis of the form 
|a;i) • • • \xn) = \xi . . . Xn), where each \xi) is |0) or |1). This enables entanglement, 
a feature necessary for the advantage of quantum computers. As an example of 
entanglement, the two-qubit state \^^) = (|00)-)-|ll))/\/2, one of the Bell states, 
can't be written as a product state \(j)i) \(j)2)- 

The linearity of quantum theory implies that a quantum computer can execute 
classical computations in superposition. For example, if the input state contains 
all possible input values |x) of a function /(x), the function can be computed 
using a unitary operation Uf as 

5^a.|x)|0)^5^a.|x)|/(x)). (1) 
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With a single call to Uf, the quantum computer produces a state that contains 
information about all the possible outputs of f{x). 

Nevertheless, quantum computation has several limitations. For example, the 



no-cloning theorem (94,63) states that an unknown quantum state cannot be 



copied perfectly. More importantly, the information of a general quantum state 
cannot be read out with a single projective measurement, because that would 
collapse a superposition into one of its components. Therefore, while the state 
in Eq. [T] contains information about all possible outputs, that information is not 
immediately accessible. Instead, a quantum algorithm has to be designed in a way 
that makes it easy to measure a global property of /, without making it necessary 
to compute all the individual f{x). Algorithms of this kind are discussed in the 
following sections. 



2.2 Approaches to quantum computing 

There are several models, or ways of formulating, quantum computation. Most 
work in quantum simulation has been done in the circuit and adiabatic models. 
While the two are known to be computationally equivalent — any computation 
that can be performed in one model can performed in the other in a comparable 
amount of time ([6,66,89) — different problems are solved more naturally in differ- 
ent models. We discuss the two models in turn, but note that other models hold 
promise for the development of future simulation algorithms, including topolog- 



ical quantum computing (69,90), one-way quantum computing (108 109), and 



quantum walks ( 65 ) . 

2.2.1 The circuit model The cornerstone of quantum computation is the 
generalization of the classical circuit model, composed of classical bits and logical 



gates, to a quantum circuit model (35,13,94). A quantum circuit is a multi-qubit 



unitary transformation U which maps a set of initial states to some final states. 
Usually, a unitary gate is decomposed into elementary gates which involve a few 
(one or two) qubits each. 
In classical computing, the nand gate is universal (95), meaning that any 



logical circuit can be constructed using nand gates only. Similarly, in quantum 
computing, there are sets of unitary operations that form universal gate sets. 
A quantum computer that can implement such a set is called universal, and 
can perform any unitary transformation U to an arbitrary accuracy. It turns 
out that the set containing all single-qubit gates in addition to any two-qubit 
entangling gate, such as cnot, is universal (63) (Fig. IT]). An entangling gate can 
be realized by any physical interaction that can generate entanglement between 
qubits. Examples of experimental implementations of quantum gates have been 
reviewed (j77p, and we will cover some of the experiments relevant to quantum 
simulation in Sec. [H 

Beside the elementary gates, an important quantum transformation is the 
quantum Fourier transform (QFT). It transforms any quantum state \ip) = 
J2x V {^) 1^) ™to its Fourier representation. 



UgPTlf) 



N-l 

E 
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V>ik)\k), 



(2) 



where 'p{k) = (l/VW) ^^=q ip{x)e^^ ' are the discrete Fourier coefficients of 
•^{x). The QFT can be efficiently implemented using a quantum circuit (94): for 
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n qubits, the number of elementary gates required is 0(71?). For comparison, the 
classical fast Fourier transform requires 0(n2") gates. We take advantage of the 



QFT in Sec. 3.2 for the simulation of quantum dynamics, and in Sec. 3.3 for the 
measurement of observables. 
2.2.2 Adiabatic quantum computation An alternative to the gate model 



is the adiabatic model of quantum computation (41 ). In this model, the quantum 
computer remains in its ground state throughout the computation. The Hamil- 
tonian H{t) of the computer is changed slowly from a simple initial Hamiltonian 
Hi to a final Hamiltonian Hj whose ground state encodes the solution to the 
computational problem. The adiabatic theorem states that if the variation of the 
Hamiltonian is sufficiently slow, the easy-to-prepare ground state of Hi will be 
transformed continuously into the ground state of Hf. It is desirable to complete 
the evolution as quickly as possible; the maximum rate of change is mostly deter- 
mined by the energy gap between the ground and first excited states during the 



evolution p6|[7| [T26 ,125). The applications of adiabatic quantum computation to 



simulation include preparing quantum states of interest and solving optimization 



problems such as protein folding (102). We discuss the details in Sees. 3.4.1 and 
|4j respectively. 

2.3 Quantum complexity theory 

To understand the computational advantages of quantum algorithms for chemical 
simulation, we discuss some aspects of computational complexity theory, which 
defines quantum speed-up unambiguously. A proper measure of the complexity of 
an algorithm is how many operations (or how much time) it takes to solve prob- 
lems of increasing size. Conventionally, a computational problem is described 
as easy or tractable if there exists an efficient algorithm for solving it, one that 
scales polynomially with input size (for an input of size n, as 0{n ) for some 
k). Otherwise, the problem is hard. This distinction is admittedly a rough one: 
for reasonable problem sizes, an "inefficient" algorithm scaling exponentially as 
0(1.0001") would be faster than an "efficient" 0{'n}^^^) algorithm. Neverthe- 
less, this convention has proven useful because, in practice, polynomially scaling 
algorithms generally outperform exponential ones. 

The class of all problema^that are easy for classical computers (classical Turing 
machines) is called P (^. The class of all problems whose answer can be verified 
in polynomial time is NP. For example, even though we don't know how to factor 
numbers efficiently, factoring is in NP because we can check the proposed answer 
efficiently by multiplication. Note that p C np because any easy problem can 
be verified easily. Whether p = np is a famously open question; however, it is 
widely believed that they are not equal, that is, that there are problems in np 
that cannot be solved easily ([T|) . The hardest among them belong to the class np- 
hard: if any NP-hard problem can be solved efficiently, then so can any problem 
in NP. 



Likewise, bqp are those problems that are easy for a quantum computer (130). 
The quantum analogue of NP is called QMA and contains those problems easy 
to check on a quantum computer. In analogy with NP-hard problems, QMA- 



hard contains the hardest problems in QMA. Shor's factoring algorithm (112) is 



^Strictly speaking, decision problems, those with a yes-or-no answer. However, other prob- 
lems can be recast as decision problems; for example, instead of asking "What is the ground-state 
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significant because it provides an example of a problem in bqp which is widely 
thought (although not proven) to be outside of P; that is, a problem believed to 
be hard on classical computers that is easy for a quantum computer. 

The relationships between the complexity classes mentioned above are illus- 
trated in Fig. [2| In the remainder of this review, we explore the advantages of 
quantum simulation over its classical counterpart, in part, by situating various 
simulation tasks in the computational classes illustrated in Fig. [2| 

3 QUANTUM SIMULATION 

Quantum simulation schemes can be divided into two broad classes. The first 
is dedicated quantum simulation, where one quantum system is engineered to 
simulate another quantum system. For example, quantum gases in optical lattices 



can be used to simulate superfluidity (12). The other, more general, approach 



is universal quantum simulation, simulating a quantum system using a universal 
quantum computeiQ Although we will focus on universal quantum simulation 
because most chemical proposals assume a universal quantum computer, we will 
mention dedicated simulators where appropriate. 

One of the main goals of quantum simulation is to determine the physical 
properties of a particular quantum system. This problem can usually be concep- 
tualized as involving three steps: 

1. Initialize the qubits in a state that can be prepared efficiently, 

2. Apply a unitary evolution to this initial statej^and 

3. Read out the desired information from the final state. 

We note at the outset that it is not possible to simulate an arbitrary unitary 
evolution on a quantum computer efficiently. An arbitrary unitary acting on a 
system of n spins has 2" x 2" free parameters, and would require an exponential 
number of elementary quantum gates to implement. However, in quantum chem- 
istry, it is usually not necessary to simulate arbitrary dynamics, since natural 
systems aren't arbitrary (j84j). Instead, the interactions between, say, molecular 
orbitals are local — featuring at most A;-body interactions — and this crucial as- 
pect of their structure can be exploited for their efficient simulation. That is, the 
Hamiltonian generating the unitary evolution is a sum H = ^^ Hi of polynomi- 
ally many terms, each of which acts on at most polynomially many degrees of 
freedom. A local Hamiltonian generates a time-evolution that can be decomposed 
into t/5t time-steps according to the Lie- Trotter formula, 

/ \ t/st 

The approximation can be improved by increasing the number of time steps or 



by using higher-order generalizations of this formula (54,18). Finally, since each 
factor e"*-'^' acts on only a sub-region of the Hilbert space and can therefore 
be efficiently simulated, so can a product of polynomially many such factors. 
Hence, the time it takes to perform the simulation scales polynomially with the 



^The terms "analog" and "digital" have also been used for dedicated and universal quantum 



simulation, respectively (24 1 

3 



Non-unitary open-system dynamics have been studied as well (111 
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simulated time t. Most methods of quantum simulation make use of the Trotter 
decomposition, and we will describe in more detail their applications in chemistry. 
We will not discuss all the available methods, for which the reader is directed to 



comprehensive reviews (27,11022). 

In the following, we describe two ways in which chemical wavefunctions can be 
encoded on a quantum computer, second- and first-quantization approaches (see 
Table IT] for a comparison). For each approach, we outline the methods of prepar- 
ing certain classes of initial states and propagating them in time. Afterward, we 
discuss the methods of measurement of observables and preparation of ground 
and thermal states, which do not depend essentially on the way the wavefunction 
is encoded. 

3.1 Second quantization 

We start by considering the purely electronic molecular problem, in which the 
Born-Oppenheimer approximation has been used to separate the electronic and 
nuclear motion. The wavefunction of the electrons can be expanded in an or- 
thonormal basis of M molecular spin-orbitals {|xi)}- Corresponding to this basis 
are the fermionic creation and annihilation operators aj and Oj. There is a very 
natural mapping between the electronic Fock space and the state of M qubits: 
Having qubit i in the state |0) (or |1)) indicates that spin-orbital i is unoccupied 
(or occupied). 

An important subtlety is that electrons in a molecule, unlike the individually 
addressable qubits, are indistinguishable. Put differently, while the operators a] 
and Oj obey the canonical fermionic anticommutation relations, {aija'A = 6ij, 
the qubit operators that change |0) to |1) and vice versa do not. This problem 
can be solved by using the Jordan- Wigner transformation to enforce the correct 



commutation relations on the quantum computer (99 118, 119, 10). The Jordan- 
Wigner transformation for this case results in the following mapping between the 
fermionic operator algebra and the qubit spin algebra: 

where <t" = |1) (0| and cr+ = |0) (1|. 

The electronic Hamiltonian in the second-quantized form is 

H = y ^ ripqapaq -|- — ^ ^ hpyrsCpagaj-cisj (5j 

pq pqrs 

where the spin-orbital indices p, q, r, s each range from 1 to M. Here, one-electron 
integrals hpq = {p\ (T -|- V/v) \q) involve the electronic kinetic energy T and the 
nuclear-electron interaction V^, and the two-electron integrals hpqrs = {pq\ Ve {rs) 
contain the electron-electron interaction term Ve- For simulation on a quantum 
computer, this Hamiltonian is recast into the spin algebra using the Jordan- 
Wigner transformation, Eq. |4j and the time-evolution it generates is implemented 
using the Trotter decomposition, Eq. ^ Note that H contains 0{M'^) terms, and 
each of these terms generates a time evolution of the form 

Tj __ -ihpqalaqSt tj __ —ihpqraa],a\aras5t /-g\ 
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Each of these operators requires 0{M) elementary quantum gates to implement 
because of the Jordan- Wigner transformation. Since there are altogether 0{M'^) 
terms that need to be implemented separately, the total cost of simulating H 
scales as 0{M^) pUj ). 

While any basis {|xi)} can be chosen to represent if , it is desirable to choose 
a basis as small as possible that adequately represents the system under study. 
Electronic structure experience provides for many good starting points, such as 
the Hartree-Fock basis or the natural orbitals (33). No matter which basis is 



chosen, a lot of the computation can be carried out on classical computers as 
pre-processing. In particular, the coefficients hpq and hpq^s can be efficiently 
pre-computed on classical computers. That way, only the more computationally 
demanding tasks are left for the quantum computer. 

Using the Hartree-Fock basis allows us to use the Hartree-Fock reference state 
as an input to the quantum computation (10). A salient feature is that such 
states are Fock states, which are easy to prepare on the quantum computer: 
some qubits are initialized to |0) and others to |1). In fact, any single-determinant 
state can be easily prepared in this way. Furthermore, it is possible to prepare 
superpositions of Fock basis states as inputs for the quantum computation. While 
an arbitrary state might be difficult to prepare, many states of interest, including 
those with only polynomially many determinant contributions, can be prepared 
efficiently (99,118,119,127). The problem of preparing an initial state that is 



close to the true molecular ground state is addressed in Sec. 3.4 



The chief advantage of the second-quantization method is that it is frugal with 
quantum resources: only one qubit per basis state is required, and the integrals 
can be pre-computed classically. For this reason, the first chemical quantum 
computation was carried out in second quantization (see Sec. pi). Nevertheless, 
there are processes, such as chemical reactions, which are difficult to describe 
in a small, fixed basis set, and for this we turn to discussing first-quantization 
methods. 



3.2 First quantization 



The ffist-quantization method, due to Zalka ( 138 , 132 , 62|, simulates particles 
governed by the Schrodinger equation on a grid in real spac^ For a single particle 
in one dimension, space is discretized into 2" points, which, when represented 
using n qubits, range from |0 . . . 0) to |1 . . . 1). The particle's wavefunction can be 



expanded in this position representation as 
to be simulated is 



H = T + V 



P 
2m. 



+ V{x), 



ar \x). The Hamiltonian 



(7) 



and the resulting unitary can be implemented using the quantum version of the 



split-operator method (^42 , 76 ) : 



U{6t) 



-iiT+V)5t 



f/QFTe-'™C// 



iVSt 



QFT'= 



(8) 



The operators e~ and e~* are diagonal in the position and momentum 
representations, respectively. A diagonal operator can be easily implemented be- 
cause it amounts to adding a phase e~*^(^)"* to each basis state |x). Furthermore, 



We also note the method of (21 1, which in our terminology is a hybrid between second- and 



first-quantization methods, ft associates a qubit to the occupation of each lattice site. 
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it is easy on a quantum computer to switch between the position and momentum 
representations of a wavefunction using the efficient quantum Fourier transform. 
Therefore, simulating a time evolution for time t involves alternately applying 
^-iVSt g^j^^ ^-iTSt -^j^j-^ ^j-^g time steps 5t chosen to be sufficiently short to secure 
a desired accuracy. Finally, the scheme can be easily generalized to many parti- 
cles in three dimensions: a system of B particles requires 2>Bn qubits, n for each 
degree of freedom. 

The first-quantization method can be applied to many problems. The earliest 
applications established that as few as 10-15 qubits would be needed for a proof- 
of-principle demonstration of single-particle dynamics (120[) (later improved to 6- 



10 ( 15 )). The method could also be used to faithfully study the chaotic dynamics 
of the kicked rotor model (jSTp. The first chemical application was the proposal 



of a method for the calculation of the thermal rate constant (|83|) (see Sec. 3.3) 



We investigated the applicability of the first-quantization method to the sim- 



ulation of chemical dynamics ( 62 ) . The simplest approach is to consider all the 
nuclei and electrons explicitly, in which case the exact non-relativistic molecular 
Hamiltonian reads 

where rij = |rj — rj| is the distance between particles i and j, which carry charges 
Qi and Qj, respectively. As before, the split-operator method can be used to sep- 
arate the unitaries that are diagonal in the position and momentum bases. Note 
that a Jordan- Wigner transformation is not required; H preserves permutational 
symmetry, meaning that if the initial state is properly (anti-)symmetrized (see 
below), it will stay so throughout the simulation. 

Since the Born-Oppenheimer approximation (BOA) has been widely used in 
quantum chemistry, it might seem extravagant to explicitly simulate all the nuclei 
and electrons. Nevertheless, the exact simulation is, in fact, faster than using the 
BOA for reactions with more than about four atoms (62) . The reason for this 
is the need to evaluate the potential l^(ri, . . . , r^) on the fly on the quantum 
computer. In the exact case, the potential is simply the pairwise Coulomb in- 
teraction; on the other hand, evaluating the complicated, many-body potential 
energy surfaces that are supplied by the BOA is a much more daunting task, even 
considering that one can use nuclear time-steps that are about a thousand times 
longer. That is, exact simulation minimizes arithmetic, which is the bottleneck of 
the quantum computation; by contrast, the bottleneck on classical computers is 
the prohibitive scaling of the Hilbert space size, which is alleviated by the BOA. 

In order to carry out simulations, it is important to prepare suitable initial 
states. Zalka's original paper (|138) contained a very general state-preparation 



scheme, later rediscovered (48,64 71) and improved (114). The scheme builds 
the state one qubit at a time by performing a rotation (dependent on the pre- 
vious qubits) that redistributes the wavefunction amplitude as desired. For 
example, Gaussian wavepackets or molecular orbitals can be constructed effi- 
ciently. We discussed how to combine such single-particle wavefunctions into 
many-particle Slater determinants, superpositions of determinants, and mixed 



states in (129). In particular, the (anti-)symmetrization algorithm of (i2j) was 
improved and used to prepare Slater determinants necessary for chemical sim- 
ulation. Furthermore, we outlined a procedure for translating states that are 
prepared in second-quantization language into first-quantized wavefunctions, and 
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vice versa. Techniques for preparing ground and thermal initial states are dis- 
cussed in Sec. 13.41 

The first-quantization approach to quantum simulation suffers from the fact 
that even the simplest simulations might require dozens of qubits and millions of 



quantum gates (62). Nevertheless, it has advantages that would make it useful 
if large quantum computers are built. Most importantly, because the Coulomb 
interaction is pairwise, simulating a system of B particles requires 0{B'^) gates, a 
significant asymptotic improvement over the second-quantized scaling of 0{M^), 
where M is the size of the basis set. 

3.3 Measuring observables 

We have discussed how to prepare and evolve quantum states on a quantum 
computer. Information about the resulting state must be extracted in the end; 
however, full characterization (quantum state tomography) generally requires 
resources that scale exponentially with the number of qubits. This is because a 
measurement projects a state into one consistent with the measurement outcome. 
Because only a limited amount of information can be extracted efficiently, one 
needs a specialized measurement scheme to extract the desired observables, such 
as dipole moments, correlation functions, etc. 

In principle, an individual measurement can be carried out in any basis. How- 
ever, since experimental measurement techniques usually address individual qubits, 
a method is needed to carry out more complicated measurements. In particular, 
in order to measure an observable A, one would like to carry out a measurement 
in its eigenbasis {|efc)}. This is achieved by the phase estimation algorithm (PEA) 
(68||3j): 



where Ck = (e/dV') and Ak are the eigenvalues of A; C—U is the unitary U = 
exp{—iAt) controlled by the ancilla qubits, which are initialized in the state 
(l/viV) ^j |t). When measuring the ancilla, the eigenvalue A^ will be measured 
with probability \ck\ and, if the eigenstates are non-degenerate, the wavefunc- 
tion will collapse to the eigenvector |efc)lj For the PEA to be efficient, it must 
be possible to simulate the pseudo-dynamics e~ efficiently. In particular, if 
we are interested in molecular energies, the observable is the Hamiltonian H, 
and we need to simulate the actual dynamics e~*^* (see Sec. ^. Note that the 
PEA is closely related to classical algorithms for preparing eigenstates by Fourier 



analysis of a propagating system (34, 122). As in classical Fourier analysis, the 



(pseudo-) dynamics must be simulated for longer in order to achieve a higher pre- 
cision in the A/^. More precisely, for a final accuracy of e, the PEA must run for 



atimeO(l/e) (94,23). 



Because quantum measurement is inherently random, repeating a measurement 
on multiple copies of the same system helps in determining expectation values 
of observables. The central limit theorem implies that measuring A^ copies of a 
state results in a precision that scales as l/\/iV (the standard quantum limit, 
SQL). For example, repeating the PEA gives an SQL estimate of the coefficients 



'Other methods for eigenvalue measurement include pairing adiabatic quantum evolution 



with Kitaev's original scheme (19 1 and applications of the Hellmann-Feynman theorem (971 
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|cfcp; these can be used to calculate the expectation value {A) = J2f^ \ck\'^Ak, also 
to the SQL. When entanglement is available, one can achieve precision scaling as 



1/N — this is the Heisenberg or quantum metrology limit (47). An algorithm for 



the expectation values of observables has been proposed that can get arbitrarily 



close to the Heisenberg limit (72). 

The first algorithm for measuring a chemical observable was Lidar and Wang's 
calculation of the thermal rate constant by simulating a reaction in first quantiza- 



tion and using the PEA to obtain the energy spectrum and the eigenstates ( 83 ) . 
These values were used to calculate the rate constant on a classical computer 
by integrating the flux-flux correlation function. We improved on this method 
with a more direct approach to the rate constant fi62j). We showed how to ef- 
ficiently obtain the product branching ratios given different reactant states — if 



the initial state is a thermal state (see Sec. 3.4.2), this gives the rate constant 
directly. Furthermore, the method was used to obtain the entire state-to-state 
scattering matrix. A method for reaction rates using a dedicated quantum sim- 
ulator where artificial molecules are experimentally manipulated, has also been 



proposed (113) 



More generally, correlation functions provide information about a system's 
transport and spectroscopic properties. On a quantum computer, the correla- 
tion function of any two observables can be estimated efficiently if their pseudo- 



dynamics can each be simulated efficiently (99 119). The method does not suffer 
from the dynamic sign problem that plagues classical Monte Carlo methods for 
sampling correlation functions. An alternative approach is the measurement of 



correlation functions using techniques of linear-response theory (124). 

Molecular properties such as the dipole moment or the static polarizability 
are also of chemical interest. They are derivatives of the molecular energy with 
respect to an external parameter, such as the electric field. We showed how to 



calculate them (61) using the PEA and the quantum gradient algorithm (59). 
The algorithm is insensitive to the dimensionality of the derivatives, an obstacle 
to classical computers. For example, the molecular gradient and Hessian can 
be computed — and used to optimize the geometry — with a number of energy 
evaluations independent of system size. 

3.4 Preparing ground states and thermal states 



In Sees. |3.1| and |3.2[ we discussed the preparation of various initial states for 
quantum simulation. We postponed discussing the preparation of ground and 
thermal states because of subtleties to which we now turn. 
3.4.1 Ground STATE PREPARATION BY PHASE ESTIMATION A large part of 
quantum chemistry is concerned with the calculation of ground-state properties 
of molecules, making it desirable to prepare such states on a quantum computer. 
In the previous section, we described how the PEA can be used to measure a 
quantum state in the eigenbasis of a Hermitian operator. This suggests a method 
for preparing a ground state: measuring in the eigenbasis of the Hamiltonian will 
project a state \ip) to the ground state \g) with probability KV'b)! • 

The problem, therefore, is to prepare a state close to the ground state, from 
which we can project the ground-state component. Choosing a random state 
iV'rand) is bound to fail, since the overlap is expected to be exponentially small 
in the number of qubits n: (^randls') ~ 2""". This means that one would have 
to repeat the PEA exponentially many times before chancing upon the ground 
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state. 

Methods of quantum chemistry can be used to improve the overlap. We studied 
the ground-state preparation of H2O and LiH in second quantization, based on 
the Hartree-Fock (HF) approximation ( |10[ ). The goal was to prepare the ground 
state of the full configuration interaction (FCI) Hamiltonian, so that its energy 
could be read out by the PEA, thus solving the electronic structure problem. 
Since these molecules were considered at equilibrium geometries, the HF guess 
was sufficient for the algorithm to estimate the ground-state energies of these 
molecules with high probability. The overlap can be further improved by choosing 
a more sophisticated approximation method such as a multi-configuration self- 



consistent field (MCSCF) wavefunction (128). 

Alternatively, the overlap can be increased using adiabatic quantum computa- 



tion (Sec. 2.2.2). We applied adiabatic state preparation (ASP) to the case of the 
hydrogen molecule H2 in the ST0-3G basis at various bond lengths (10). As the 



bond length increases, the HF state has decreasing overlap with the exact state, 
reaching 0.5 at large separations. ASP works by preparing the ground state of the 
HF Hamiltonian and then slowly changing to the FCI Hamiltonian. The speed 
of the variation of the Hamiltonian is limited by the energy gap between the 
ground state and the first excited state. In the case of H2, this method allowed 
the preparation of the FCI ground state with a high fidelity. 

Procedures similar to ASP have been proposed to study low-energy states of 



some toy models in physics (97) and superconductivity (134). It is also possible 
to encode a thermal state into the ground state of a Hamiltonian (|5| |116 ) , offering 
a way to prepare a thermal state, a problem further discussed in the next section. 
3.4.2 Thermal state preparation While not often a subject of quantum- 
chemical calculations, the thermal states are of significance because they can be 
used to solve many problems, ranging from statistical mechanics to the calcu- 
lation of thermal rate constants. Classical algorithms typically rely on Markov 
chain Monte Carlo (MCMC) methods, which sample from the Gibbs density ma- 
trix, p = e~^^ /Z, where Z is the partition function. The challenge is that it is 
generally impossible to sample from the eigenstates le^) of a certain Hamiltonian 
H if they are not pre-determined (which is often more challenging) . 

With a quantum computer, assuming the PEA can be efficiently implemented, 
we can prepare the thermal state of any classical or quantum Hamiltonian from a 



Markov chain constructed by repeating a completely positive map (124). A limi- 



tation of that approach is that the Metropolis step can make too many transitions 
between states of very different energies, sometimes leading to a slow convergence 
rate of the resulting Markov chain. This issue was addressed by building up the 



Markov chain by applying random local unitary operations (123). The result 



ing operation is a Metropolis- type sampling for quantum states; although the 
underlying Markov chain is classical in nature, performing it on a quantum com- 
puter provides the benefit of being able to use the PEA without explicitly solving 
the eigenvalue equations. However, quantum computers can implement Markov 
chains corresponding to thermal states of classical Hamiltonians with a quadratic 



speed-up (121, 117 133, 106, 107 26). 



Zalka's state preparation algorithm (see Sec. 3.2) is applicable to preparing the 



coherent encoding of thermal states (CETS) 1-0 



th/ 



\Ai.) = J2\f^^^^\^k)\ek) 



:n) 
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which is equivalent to the Gibbs density matrix, pth = Z^fc e^^^^^/Z |efc) (efc|, if 
one register is traced out. If the eigenstates and eigenvalues are known, it is 



possible to construct the GETS directly (129). On the other hand, combining 
ideas from belief-propagation ( J87| and quantum amplitude amplification (l63l), 
we were able to construct the GETS of classical Hamiltonians with a quadratic 



quantum speedup (137) 



Lastly, a thermal state can be prepared by modeling the physical interaction 



with a heat bath (138, 124). However, the computational cost of these methods 
is not well understood. 
3.4.3 QMA-HARDNESS AND FUTURE PROSPECTS Unfortunately, the pro- 
cedures for ground- and thermal-state preparation outlined above are not fully 
scalable to larger systems. A simple way to see this is to imagine a system 
composed of N identical, non-interacting molecules. Even if one molecule can 
be prepared with a ground-state overlap of 1 — e by any method, the fidelity of 



the A'^-molecule state will be exponentially small, (1 — e) (73). ASP would fail 
when the energy gap got so small that the Hamiltonian would have to be varied 
exponentially slowly. 

More generally, there are broad classes of Hamiltonians for which finding the 
ground state energy (and therefore also a thermal state) is known to be QMA- 



hard, that is, most likely hard even on a quantum computer (see Sec. 2.3) (70 



4,66 98,111). Nevertheless, the scaling of the ground- and thermal-state energy 
problems for chemical systems on a quantum computer is an open question. It 
is possible that algorithms can be found that are not efficient for all QMA-hard 
Hamiltonians, but nevertheless succeed for chemical problems. 

4 OPTIMIZATION WITH ADIABATIC QUANTUM SIMU- 
LATION 

We describe the use of quantum computers to solve classical optimization prob- 
lems related to chemistry and biology. This class of problems plays an important 
role in fields such as drug design, molecular recognition, geometry optimization, 



protein folding, etc. (53,44). 

Of all the models of quantum computation, adiabatic quantum computation 
(AQG) is perhaps the best suited for dealing with discrete optimization prob- 



lems. As explained in Sec. 2.2.2 the essential idea behind AQG is to encode the 
solution to a computational problem in a (final) Hamiltonian ground state which 
is prepared adiabatically. 

Although final Hamiltonians have been proposed for various problems related 
to computer science ( 4Tp6| 93 .92,29]), only recently we derived constructions ( 102 ) 



for problems of chemical interest such as the lattice heteropolymer problem ( 100 
88 , 75 ) , an NP-hard problem ( [52] ) . It can be used as a model of protein folding 



(37), one of the cornerstones of biophysics. Note that the quantum-computational 
implementation of the protein folding problem does not assume that the pro- 
tein is treated quantum mechanically. Instead, the quantum computer is being 
used as a tool to solve the classical optimization problem. In the lattice folding 
problem, the sequence of amino acids is coarse-grained to a sequence of beads 
(amino acids) connected by strings (peptide bond). This chain of beads occupies 
points on a two- or three-dimensional lattice; a valid configuration (fold) is a 
self-avoiding walk on the lattice and its energy is determined by the interaction 
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energies among amino acids that are non-bonded nearest neighbors in the lattice. 



The hydrophobic-polar (HP) model (80) is the simplest realization of this prob- 



lem. The amino acids are broken into two groups, hydrophobic (H) and polar (P). 
Whenever two non-bonded H amino acids are nearest neighbors in the lattice, 
the free-energy of the protein is reduced by one unit of energy, Ehh = — 1- The 
remaining interactions do not contribute to the free energy Ehp = Epp = 0. The 
lattice folding problem consists in finding one of more folds that minimize the 
free energy of the protein. By the thermodynamic hypothesis ( |39[ ), such fold(s) 
correspond to the conformation of the native conformation(s) of the protein. 

The theory behind the quantum-computational implementation of lattice fold- 
ing is guided by the proposed quantum adiabatic platform on superconducting 



qubits (60). This scheme is designed to find solutions to the problem. 



E{s) = - ^ hiSi + ^ JijSiSj, (12) 

i j>i 

where |/ij| < 1, \Jij\ < 1, and Si = ±1. Given a set of {hi} and the interaction 
matrix {Jij}, the goal is to find the assignment s* = s^s"^ ■ ■ ■ s*j^ , that minimizes 
E{s). 

The time-dependent Hamiltonian is chosen to be, 

H{r) = A{r)H, + B{T)Hf, r = t/tmn, (13) 

where Hi has a simple-to-prepare ground state and Hf = — J2i ^«'^2+Z^}>j Jij'^l'^z, 
with (T* denotes the Pauli matrix acting on the ith qubit, and tmn is the running 
time. The time-dependent functions A{t) and B{t) are such that A(0) ^ ^(0) 
and A{1) <^ -B(l). Therefore, at the beginning (end) of the simulation, the 
ground state corresponds to the ground state of Hi {Hf). Note that, as desired, 
|s*) = |sj, ^2, • • • , s^) is the ground state of Hf. Measurement of this final state 
provides the solution to our problem. 

The theoretical challenge is to map the lattice folding free energy function into 



the form of Eq. 12 |102, 104). In two dimensions, we use two binary variables 
determining the direction of each bond between two amino acids (beads). If a 
particular bond points upwards, we write "11"; if it points downwards, leftwards 
or rightwards, we write "00", "10", or "01", respectively. For an A^-amino-acid 
protein, we need two binary variables for each of the A^ — 1 bonds. Fixing the 
direction of the first bond reduces the number of variables to £ = 2(A^ — 2) binary 
variables. Any possible A^-bead fold can be represented by the string of binary 
variables of the form 01gig2 • • • Qe-iQe, where we set the direction of the first bond 
to be right ("01"). 



As an example, the free energy function (104) associated with the folding of a 



four-amino-acid peptide assisted by a "chaperone" protein (see Fig. [3|) is 

E{q) = 4 - 3gi -F 4^2 - 4gig2 - qs + qiQs - 2g2g3 + 40^4 - 2gig4 , . 

- 8^234 + 5^15294 - 2^394 + 5q2q3q4 - qiq^qmi- 

By substituting values for the four binary variables defining the directions of 
the second (^1^2) and third ((73(74) bonds, we can verify that the 16 assignments 
provide the desired energy spectrum (Fig. ^. Eq. 14 is not in the form of Eq. [12} 



We converted this energy function fr om i ts quartic form to a quadratic form, using 
two extra ancilla binary variables (104). After the substitution qi = 2(1 — Sj), 
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the free energy function now resembles that of Eq. 12 
reahzation is described in Sec. [5j 



An early experimental 



Note that solving the HP model is NP-hard (16,32,52). AQC is equivalent to 



the circuit model, so it is unlikely able to solve NP-hard problems in polynomial 
time (see Sec. 2.3). Real world problems (and the instances defining biologi- 
cally relevant proteins) are not necessarily structureless. Taking advantage of the 
structure or information about a particular problem instance is one of the ideas 
behind new algorithmic strategies (8,40,105). An example is to introduce heuris- 



tic strategies for AQC by initializing the calculation with a educated guess ( 105 ) 



5 EXPERIMENTAL PROGRESS 



Experimental quantum simulation has rapidly progressed (24,22) since the early 



simulation of quantum oscillators using nuclear magnetic resonance (NMR) (115). 
Here we review chemical applications of available quantum-computational de- 
vices. 

Quantum optics On an optical quantum computer, various degrees of free- 
dom of single photons, such as polarization or path, are used to encode quantum 
information ( 96 , 74 ) . This architecture was used for the first quantum simula- 



tion of a molecular system, a minimal basis model of the hydrogen molecule H2 
(79). Qubits were encoded in photon polarization, while two-qubit gates were 



implemented probabilistically using linear-optical elements and projective mea- 
surement. The minimal-basis description of H2 used two spin-orbitals per atom. 
Since the FCI Hamiltonian is block-diagonal with 2x2 blocks, two qubits suf- 
ficed for the experiment: one for storing the system wavefunction, and one for 
the readout of the PEA. The PEA was implemented iteratively, extracting one 
bit of the value of the energy at a time. Twenty bits of the energy were obtained, 
and the answer was exact within the basis set. Fig. |4] describes the experiment 
and the potential energy surfaces that were obtained. 
Nuclear magnetic resonance Nuclear spins can serve as qubits, being 



addressed and read-out using an NMR spectrometer ( 14 ) . The first experimental 



quantum simulation, of a harmonic oscillator, was performed using NMR (115). 



The platform has since been used to simulate a number of model systems (91,23 
135 101 ), leading up to the recent simulation of H2 (f38J). The H2 experiment used 
^''C-labeled chloroform, in which the carbon and hydrogen nuclear spins form two 
qubits. The experiment achieved 45 bits of precision (15 iterations of PEA, 3 bits 



per iteration) in the ground state energy. Adiabatic state preparation (Sec. 3.4) 
was implemented for various bond distances. 

Superconducting systems The circulating current (clockwise or counter- 
clockwise) flowing in a micron-sized loop of a superconductor can be used as a 
qubit ( 136[[82 ). Examples of applications based on superconducting qubits include 
the tailor-made generation of harmonic oscillator states ( 55 ) and the implemen- 
tation of the Deutsch- Jozsa and Grover quantum algorithms ( 36 ) . Recently, the 



free energy function discussed in Sec. |4] for the four-amino-acid peptide assisted 
by a chaperone protein (see Fig. ^ has been experimentally realized ( 103| ). A 
microprocessor consisting of an array of coupled superconductor qubits has been 
used to implement the time-dependent Hamiltonian in Eq. 13, with Hi oc X^i "^x 



as the initial Hamiltonian ( 60 , 58 , 49 ) . The quantum hardware operating at a 



temperature of 20 mK found the correct solution with a probability of 



Simulating Chemistry using Quantum Computers 15 

Characterization of this device is currently underway ([5l}[l7|[50|[78| . 
Trapped ions Qubits can also be encoded in the electronic states of cold 



trapped ions, offering one of the most controllable systems available today (30 



57 20). This platform has already produced sophisticated simulations of physical 



systems (45,67,46), but chemical applications are still to come. 



6 CONCLUSIONS 

We have outlined how a quantum computer could be employed for the simulation 
of chemical systems and their properties, including correlation functions and re- 
action rates. A method for lattice protein folding was also discussed. Although 
we focused on the adiabatic and circuit models, these are not the only universal 
models of quantum computation and it may be possible to make further algo- 



rithmic progress with models such as topological quantum computing ( 69 90 ) , 



one-way quantum computing (108,109), and quantum walks (65,28,85). 

We reported on the first experiments relevant to chemistry and we expect 
more to come in the near future. With recent technological advances, there are 
many prospects for the future of quantum simulation. However, as more qubits 
are added to experiments, more effort will be needed to control decoherence, 
since error correction procedures ( j94p might not be sufficient in practice due to 



the spatial and temporal overheads required (31). Instead, it may be possible to 
build resilient quantum simulators or to incorporate the noise into the simulation. 

Although practical quantum computers are not available yet, quantum informa- 
tion theory has already influenced the development of new methods for quantum 
chemistry. For instance, density matrix renormalization group has been extended 
using quantum information, and its applications to chemistry have been vigor- 
ously pursued (|25j) . By studying the simulation of chemical systems on quantum 
computers, we can also expect new insights into the complexity of computing 
their properties classically. 

In analogy to classical electronics, one could say that, as of 2010, the im- 
plementation of quantum information processors is in the vacuum-tube era. A 
development parallel to that of the transistor would allow for rapid progress in the 
capabilities of quantum information processors. These larger devices would allow 
for routine execution of exact, non-adiabatic dynamics simulations, as well as full- 
configuration interaction calculations of molecular systems that are intractable 
with current classical computing technology. 

Summary Points 

• A universal quantum computer can simulate chemical systems more effi- 
ciently (in some cases exponentially so) than a classical computer. 

• Preparing the ground state of an arbitrary Hamiltonian is a hard (qma- 
complete) problem. However, the ground state of certain chemical Hamil- 
tonians can be found efficiently using quantum algorithms. 

• Simulation of quantum dynamics of physical systems is in general efficient 
with a quantum computer. 

• Properties of quantum states can be obtained by various measurement 
methods. 
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• Classical optimization problems, such as lattice protein folding, can be stud- 
ied by means of the adiabatic quantum-computational model. 

• Quantum simulation for chemistry has been experimentally realized in quan- 
tum optics, nuclear magnetic resonance, and superconducting devices. 

Future Issues 

• Developing quantum simulation methods based on alternative models of 
quantum computation is an open research direction. 

• Dedicated quantum simulators built so far are mostly for simulating con- 
densed matter systems. It is desirable to make experimental progress on 
simulating chemical systems. 

• Decoherence is currently the major obstacle for scaling up the current ex- 
perimental setups. Progress in theoretical and experimental work is needed 
to overcome decoherence. 

• We have not covered methods of quantum error correction, which will be 
important for large scale simulations. 
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Figure 1: Qubit, elementary gates and phase estimation algorithm: (a) 

The quantum state of a qubit can be represented in a Bloch sphere, (b) The action 
of the Hadaniard gate H on a qubit is shown in (a). The CNOT gate, together 
with single qubit gates, form a universal gate set. The quantum circuit for the 
phase estimation algorithm PEA is shown on the right. Here U \(j)k) = e^'^*'^'-' \(pk) 
and QFT is the quantum Fourier transform. The eigenvalues in this case have 
4-digit accuracy. 




Figure 2: Computational complexity classes. Shown are the conjectured 
relationships between the computational complexity classes discussed in this re- 
view. Simulating the time-evolution of chemical systems (denoted by the star) is 
in BQP but widely believed to be outside of P (assuming a constant error and sim- 
ulated time). That is, it is easy on quantum computers, but probably hard — even 
in principle — on conventional ones. 
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Advantages 



• Compact wavefunction repre- 
sentation (requires fewer qubits) 

• Takes advantage of classical 
electronic-structure theory to 
improve performance 

• Already experimentally imple- 
mented 



Better asymptotic scaling (re- 
quires fewer gates) 

Treats dynamics better 

Can be used for computing reac- 
tion rates or state-to-state tran- 
sition amplitudes 



Table 1: Comparison of second- and first-quantization approaclies to 
quantum simulation. 
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Figure 3: Lattice folding using quantum adiabatic hardware: In vivo 
folding of proteins involves the assistance of molecular chaperone proteins whose 
main function is to assist the folding of the newly synthesized polypeptide, (a) 
In vacuo, the four-amino-acid peptide could fold clockwise or counterclockwise, 
depending if the third amino acid moves downwards or upwards, respectively, 
(b) The chaperone molecule (represented as the pink region) obstructs the third 
amino acid from moving downward, making the counterclockwise folding the only 



global minima in the energy landscape, (c) Energy landscape (Eq. 14): Each 



overlap of the amino acids with the chaperone raises the free energy by four 
units, while overlap among amino acids in the chain raises the free energy by 
two units. The four binary variables encode the direction of the third and fourth 
bond of the peptide (see text). The quantum adiabatic hardware found the right 
solution 78 % of the times (|^3J). 
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Figure 4: Experimental simulation of the H2 molecule on a linear-optical 



quantum computer (79). a) Two-qubit iterative version of the phase esti- 
mation algorithm for evaluating molecular energies, b) Decomposition of the 
algorithm into gates, c) The layout of the optical elements used to implement 
the quantum gates on photonic polarization qubits. d) The computed potential 
energy surfaces of the H2 molecule in a minimal basis set. The results are the 
exact (in the basis) full configuration interaction energies, to 20 bits of precision. 



